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Abstract 

Simple models of a classical spacetime foam are considered, which consist of identical static de- 
fects embedded in Minkowski spacetime. Plane-wave solutions of the vacuum Maxwell equations 
with appropriate boundary conditions at the defect surfaces are obtained in the long-wavelength 
limit. The corresponding dispersion relations u> 2 = uj 2 {k) are calculated, in particular, the coeffi- 
cients of the quadratic and quartic terms in k. Astronomical observations of gamma-ray bursts and 
ultra-high-energy cosmic rays then place bounds on the coefficients of the dispersion relations and, 
thereby, on particular combinations of the fundamental length scales of the static spacetime-foam 
models considered. Spacetime foam models with a single length scale are excluded, even models 
with a length scale close to the Planck length (as long as a classical spacetime remains relevant). 
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I. INTRODUCTION 



Whether or not space remains smooth down to smaller and smaller distances is an open 
question. Conservatively, one can say that the typical length scale of any fine-scale structure 
of space must be less than 10 -18 m = 10~ 3 fm, which corresponds to the inverse of a center- 
of-mass energy of 200 GeV in a particle-collider experiment. Astrophysics provides us, of 
course, with much higher energies, but not with controllable experiments. Still, astrophysics 
may supply valuable information as long as the relevant physics is well understood. 

In this article, we discuss astrophysical bounds solely based on solutions of the Maxwell 
(and Dirac) equations. These solutions hold for a particular type of classical spacetime 
with nontrivial small-scale structure. Specifically, we consider a static (time-independent) 
fine-scale structure of space, which is modeled by a homogeneous and isotropic distribution 
of identical static "defects" embedded in Minkowski spacetime. With appropriate boundary 
conditions at the defect surfaces, plane-wave solutions of the vacuum Maxwell equations are 
obtained in the long- wavelength limit. That is, the wavelength A must be much larger than 
max(6, /), with b the typical size of the individual defect and I the mean separation between 
the different defects. An (imperfect) analogy would be sound propagation in a block of ice 
with frozen-in bubbles of air. 

Generalizing the terminology of Wheeler and Hawking Aaaaaa , we call any classical 
spacetime with nontrivial small-scale structure (resembling bubbly ice, Swiss cheese, or 
whatever) a "classical spacetime foam." The plane- wave Maxwell solutions from our classical 
spacetime-foam models, then, have a modified dispersion relation (angular frequency squared 
as a function of the wave number k = \k\ = 2tt/X): 



[«,« , y pc r] = M c 2 ,.2 + M ( )2 ^ ^ + ^ (1.1) 

where c is the characteristic velocity of the Minkowski line element (ds 2 = c 2 dt 2 — \dx\ 2 ), 
a^i and dimensionless coefficients depending on the fundamental length scales of 

the model (one length scale being b), and r labels different kinds of models. 

For simplicity, we consider only three types of static defects (or "weaving errors" in the 
fabric of space): 

(i) a nearly pointlike defect with the interior of a ball removed from R 3 and antipodal 
points on its boundary identified; 

(ii) a nearly pointlike defect with the interior of a ball removed from R 3 and boundary 
points reflected in an equatorial plane identified; 

(iii) a wormhole-like defect with two balls removed from R 3 and glued together on their 
boundaries; cf. Refs. 

Further details will be given in Sec. Ill Al As mentioned above, the particular spacetime 
models considered consist of a frozen gas of identical defects (types r = 1, 2, 3) distributed 
homogeneously and isotropically over Euclidean space R 3 . We emphasize that these classical 
models are not intended to describe in any detail a possible spacetime-foam structure (which 
is, most likely, essentially quantum-mechanical in nature) but are meant to provide simple 
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and clean backgrounds for explicit calculations of potential nonstandard propagation effects 
of electromagnetic waves. 

The type of Maxwell solution found here is reminiscent of the solution from the so-called 
"Bethe holes" for waveguides 0] • In both cases, the standard Maxwell plane wave is modified 
by the radiation from fictitious multipoles located in the holes or defects. But there is a 
crucial difference: for Bethe, the holes are in a material conductor, whereas for us, the 
defects are "holes" in space itself. 

Returning to our spacetime-foam models, we also calculate the modified dispersion rela- 
tion of a free Dirac particle, for definiteness taken to be a proton (mass m p ): 



fj h- 2 c 4 m 2 p + a [ ; ] 2 c 2 k 2 + a [ ;\ ( &M ) 2 c 2 k' + . . . , (1.2) 

with reduced Planck constant h = h/2n and dimensionless coefficients a p o, a P 2, and a p ^ . As 
might be expected, the response of Dirac and Maxwell plane waves to the same spacetime- 
foam model turns out to be quite different, with unequal quadratic coefficients a P 2 and 
a 7 2, for example. The different proton and photon velocities then allow for Cherenkov-type 
process es la , M- But, also in the pure photon sector, there can be interesting time-dispersion 
effects [10( as long as the quartic coefficient a 74 of the photon dispersion relation (II. ip is 
nonvanishing. 

In fact, with the model dispersion relations in place, we may use astronomical observations 
to put bounds on the various coefficients a 2 and a^, and, hence, on particular combinations 
of the model length scales (e.g., average defect size b and separation 7). Specifically, the 
absence of time dispersion in an observed TeV flare from an active galactic nucleus bounds 
|a 7 4| and the absence of Cherenkov-like effects in ultra-high-energy cosmic rays bounds 
) and a 7 4. In other words, astrophysics not only explores the largest structures of 
space (up to the size of the visible universe at approximately 10 10 lyr pa 10 26 m) but also the 
smallest structures (down to 1CT 26 m or less, as will be shown later on). 

The outline of the remainder of this article is as follows. In Sec. HJJ we discuss the cal- 
culation of the effective photon dispersion relation from the simplest type of foam model, 
with static r = 1 defects. The calculations for isotropic r = 2 and r = 3 models are similar 
and are not discussed in detail (App. [A] gives additional results for anisotropic defect distri- 
butions). Some indications are, however, given for the calculation of the proton dispersion 
relation from model r = 1 with details relegated to App. [Bl The main focus of Sec. [Til 
and the two appendices is on modified dispersion relations but in Sec. Ill CI we also discuss 
the Rayleigh-like scattering of an incoming electromagnetic wave by the model defects. In 
Sec. IHfl we summarize the different dispersion relations calculated and put the results in 
a general form. In Sec. IIVI this general photon dispersion relation is confronted to the as- 
tronomical observations and bounds on the effective length scales are obtained. In Sec. |V] 
we draw an important conclusion for the classical small-scale structure of space and present 
some speculations on a hypothetical quantum spacetime foam. 
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II. CALCULATION 



A. Defect types 

The present article considers three types of static defects obtained by surgery on the 
Euclidean 3-space R 3 . The discussion is simplified by initially choosing the origin of the 
Cartesian coordinates x = (x 1 , x 2 , x 3 ) = (x,y,z) of R 3 to coincide with the "center" of 
the defect. The corresponding Minkowski spacetime IxK 3 has standard metric (rj^) = 
diag(l, —1, —1, —1) for coordinates x M = (x°,x) = (ct,x m ) with index \i = 0, 1, 2, 3. 

The first type of defect (label r = 1) is obtained by removing the interior of a ball from 
R 3 and identifying antipodal points on its boundary. Denote this ball, its boundary sphere, 



and point reflection by 

B b = {iGl 3 : \x\ < b] , (2.1a) 

S b = [x E R 3 : |f| = b} , (2.1b) 

P(x) = -x. (2.1c) 
Then, the 3-space with a single defect centered at the origin x = is given by 

M[ r=11 = m£= 1] = {x e R 3 : \x | > b A (S b 3 x S P{x) e S b )} , (2.2) 



[r=ll 

where = denotes pointwise identification and M b is a shorthand notation. The 3-space 
(12. 2p has no boundary because of the S b identifications (cf. Figure [T]) and, away from the 
defect at |x| = b, is certainly a manifold (hence, the suggestive notation M). The resulting 
spacetime is M = R x M^ -1 ' . 

The corresponding classical spacetime-foam model is obtained from a superposition of 
r = 1 defects with a homogeneous distribution. The number density of defects is denoted 
n = l~ 3 and only the case of a very rarefied gas of defects is considered (b <S I), so that there 
is no overlap of defects. Clearly, there is a preferred reference frame for which the defects 
are static. Such a preferred frame, in the context of cosmology, may or may not be related 
to the preferred frame of the isotropic cosmic microwave background. 

In more detail, the construction is as follows. The 3-space with N > 1 identical defects 
is given by 

Mfr 1] , , . = M [ I= 1] n M [ I= 1] n . . . n M [ I =1 l , (2.3) 

where M^T - ^ is the single-defect 3-space (12.21) with the center of the sphere moved from 
x = to x = x n . The minimum distance between the different centers x n of 3-space (12. 3p 
is assumed to be larger than 2b. The final spacetime foam model results from taking the 
Cartesian product of R with the iV — > oo "limit" of (12.31) . 

-^iiSLion,, = R x (^M^^,) , (2.4) 

where one needs to give the statistical distribution of the centers x n . As mentioned above, 
we choose the simplest possible distribution, homogeneous, and the quantity to specify is 
the number density n of defects. 
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The second type of defect (r = 2) follows by the same construction, except that the 
identified points of the sphere Sb are obtained by reflection in an equatorial plane with unit 
normal vector a. For a point x on the sphere S&, the reflected point is denoted i2g(x). [With 
only one defect present, global Cartesian coordinates can be chosen so that a points in the 
3-direction and ( 00 j 00 >j 00 ) G Sb is to be identified with ( OC >j 00 j oc 3 ) G Sb] The space with 
a single r = 2 defect centered at the origin x = (not indicated by our shorthand notation) 
is then given by (cf. Figure [2]) 



Ml"" 21 = {fel 3 : \x | > b A (S b 3 x = Ra(x) eS b )}, (2.5) 
and spacetime is R x M^ 2 '. However, the defect embedded in (12.51) is not a manifold but 



an orbifold i.e., a coset space M/G, for manifold M and discrete symmetry group G. 
The 3-space (12.51) has, in fact, singular points corresponding to the fixed points of Rs, which 
lie on the great circle of Sb in the equatorial plane with normal vector a. But away from 
these singular points, the 3-space is a genuine manifold and we simply use the (slightly 
misleading) notation M in (12.51) . The corresponding classical spacetime- foam model results 
from a homogeneous and isotropic (randomly oriented) distribution of r = 2 defects. 

The third type of defect (r = 3) is obtained by a somewhat more extensive surgery 0,0]. 
Now, the interiors of two identical balls are removed from M 3 . These balls, denoted Bb and 
B' b , have their centers separated by a distance d > 2b. The two boundary spheres Sb and 
S' b are then pointwise identified by reflection in the central plane. This reflection is again 
denoted Ra- [With ball centers at x = (±d/2, 0, 0), the reflection plane is given by x 1 = 0] 
The space manifold with a single r = 3 defect centered at the origin x = is now (cf. Figure 
ED 



Mg-J = {f6l 3 : \x- (d/2) a\>bA \x+ (d/2) a\>bA(S b 3x ^ R u {x) G S' b ) } , (2.6) 

and the spacetime manifold is R x Mg . Using standard wormhole terminology BB a a, 
0] , the static r = 3 defect has two "wormhole mouths" of diameter 2b, with corresponding 
points on the wormhole mouths separated by a vanishing distance through the "wormhole 
throat" and by a "long distance" D G [d — 2b, d + nb] via the ambient Euclidean space. 
Again, the corresponding classical spacetime-foam model results from a homogeneous and 
isotropic distribution of r = 3 defects. 

For later use, we also define a r = 3' defect with distance d set to the value 46, where the 
factor 4 has been chosen arbitrarily. An individual r = 3' defect then has only one length 
scale, b, which simplifies some of the discussion later on. The relevant space manifold is thus 

M^ 3 ' 1 - m%=i , (2.7) 

in terms of the r = 3 manifold defined by (12.61) . 

Let us end this subsection with two parenthetical remarks, one mathematical and one 
physical. First, the r = 1 and r = 3 spaces are multiply connected (i.e., have noncontractible 
loops) but not the r = 2 space. Second, the classical spacetimes considered in this article do 
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not solve the vacuum Einstein equations but appear to require some exotic form of energy 
located at the defects; see, e.g., Part III of Ref. p for further discussion. In one way or 
another, the fine-scale structure of spacetime may very well be related to the "dark energy" 
of cosmology; see, e.g. Refs. 12|, |_13(] and references therein. 




FIG. 1: Three-space (|2.2p from a single spherical defect (type r = 1, radius b) embedded in M 3 , 
with its "interior" removed and antipodal points identified (as indicated by the pairs of open and 
filled circles). The corresponding classical spacetime- foam model has a homogeneous distribution 
of static r = 1 defects embedded in Minkowski spacetime. 




FIG. 2: Three-space (12. 5p from a single spherical defect (type r = 2, radius b) embedded in M 3 , 
with its "interior" removed and points identified by reflection in the equatorial plane with normal 
vector a. The corresponding classical spacetime-foam model has a homogeneous and isotropic 
distribution of static r = 2 defects embedded in Minkowski spacetime. 




FIG. 3: Three-space (12.6P from a single wormhole-like defect (type r = 3, two spheres with radii 
b and distance d between the centers) embedded in M 3 , with the "interiors" of the two spheres 
removed and their points identified by reflection in the central plane with normal vector a. The 
corresponding classical spacetime-foam model has a homogeneous and isotropic distribution of 
static r = 3 defects embedded in Minkowski spacetime. 



6 



B. Photon dispersion relation 



The task, now, is to determine the electromagnetic-wave properties for the three types of 
classical spacetime-foam models considered. The r = 1 case will be discussed in some detail 
but not the other cases (r = 2, 3), for which only results will be given. 

The calculation is relatively straightforward and consists of three steps. First, recall the 
vacuum Maxwell equations in Gaussian units 



V-£ = 0, cVxE + d t B = 0, (2.8a) 

V-£? = 0, cVxB-d t E = 0, (2.8b) 

and the standard plane-wave solution over Minkowski spacetime, 

E (x,t) = E x exp(i kz — icu y t) , (2.9a) 

B (x,t) = E y exp(i kz — ia> 7 t) , (2.9b) 

with amplitude E and dispersion relation u 2 = c 2 k 2 . This particular solution corresponds 
to a linearly polarized plane wave propagating in the z = x 3 direction (x and y are unit 
vectors pointing in the x 1 and x 2 directions, respectively). The pertinent observation, now, 
is that the electromagnetic fields (12. 9p also provide a valid solution of the Maxwell equations 
between the "holes" of the classical spacetime-foam models considered, for example, model 
(E3D for r = 1 defects. 

Second, add appropriate vacuum solutions (E\,B\), so that the total electric and mag- 
netic fields, E = Eq + E\ and B = B + B\, satisfy the boundary conditions from a single 
defect. The specific boundary conditions for the electric field at the defect surface follow by 
considering the allowed motions of an electrically charged test particle and similarly for the 
boundary conditions of the magnetic field. Geometrically, the electromagnetic-field bound- 
ary conditions trace back to the proper identification of the defect surface points and their 
tangential spaces. 

Third, sum over the contributions (Ej,Bj) of the different defects (J = 1, 2, 3, . . .) in 
the model spacetime foam and obtain the effective dielectric and magnetic permeabilities, e 
and /i, which may be wavelength dependent. The dispersion relation for the isotropic case 
is then given by 

u 2 (k)=c 2 k 2 /(e(k)^k)), (2.10) 

and we refer the reader to the textbooks for further discussion (see, e.g., Sec. II 32 3 of 
Ref. (3 and Sees. 7.5(a) and 9.5(d) of Ref. [3). 



The specifics of the second and third step of the calculation for r = 1 defects are as 
follows. In step 2, the motion of a test particle under influence of the Lorentz force (see, for 
example, the points marked in Fig. [T] for tangential motion) gives the following boundary 
conditions at the surface S& of 3-space (12. 2ft : 

- - [t=1] 

n(x) ■ E(x,t) = —n(—x) ■ E(—x,t) , n(x) x E(x,t) = +n(— x) x E(—x,t) , (2.11a) 

\x\=b 
[r=l] 

n(x) ■ B(x,t) = +n(—x) ■ B(—x,t) , n(x) x B(x,t) = — n(— x) x B(—x,t) , (2.11b) 

151=6 
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where n(x) is the outward unit normal vector of the surface at point x G Sb C R 3 . The 
boundary conditions (12.111) also ensure an equal Poynting vector c(E x B)/Att at antipodal 
points, as might be expected for an energy- flux density passing through the defect. As 
mentioned above, these boundary conditions trace back to the antipodal identification of 
the points on the surface Sb and the identification of the respective tangential spaces. 

Constant fields E oc E x and B oc E y, corresponding to the unperturbed fields (12. 9p 
over distances of order b <C A, do not satisfy the defect boundary conditions (12.111) and 
need to be corrected. As discussed in the Introduction, the correction fields Ei and Bi 
correspond to multipole fields from "mirror charges" located inside the defect. The leading 
contributions of a r = 1 defect come from fictitious electric and magnetic dipoles at the 
defect center, each aligned with their respective initial fields (12.91) and both with a strength 
proportional to b 3 E Q . For r = 1 and kb <C 1, the electric field E = Eq + E% turns out to 
be normal to the surface Sb and the magnetic field B = Bo + .Bitangential, just as for a 
perfectly conducting sphere (see, e.g., Sees. 13.1 and 13.9 of Ref. [14j]). 

In step 3 of the calculation for r = 1 defects, the effective dielectric and magnetic per- 
meabilities (Gaussian units) are found to be given by 

e [r=1] ~1 + 47rnb 3 [j (kb) + j 2 (kb)] , (2.12a) 
fj[r=n „ i _ 2nnb 3 [ 3o (kb) + j 2 (kb)} , (2- 12b) 

where n = l/l 3 is the number density of defects (mean separation I) and j p (z) is the spherical 
Bessel function of order p, for example, jo(z) = (sin z)/z. The similarity signs in Eqs. (12.12aj) 



and (12.12bl) indicate that only the i = 1 multipoles have been taken into account 171 ] . With 
(I2.10p . the dispersion relation is then 



^ {k) ) ~ (l + 4iinb 3 ij (kb)+j 2 (kb)})(l-2nnb 3 lMkb)+j 2 (kb)}) ' (2 ' 13) 

which holds for kb <C kl 1. A Taylor expansion in nb 3 and kb gives 

(c4 r=1 l (A;)) 2 ~ (1 - 2vr nb 3 ) c 2 k 2 + (vr/5) nb 5 c 2 k A + . . . . (2. 14) 

As mentioned already, the dispersion relation (I2.13P has only been derived for sufficiently 
small values of k. However, taking this expression (12.131) at face value, we note that the 



corresponding front velocity (ff ron t = lim^oo v phase) would be precisely c. See, e.g., Ref. [18 
for the relevance of the front velocity to the issues of signal propagation and causality. 
For the r = 2 spacetime-foam model, the dispersion relation is found to be given by 

J " 2,{k) Y ~ i+ ^ £L) ■ < 2 ' i5 > 

with Taylor expansion 

(^ r=2l (£;)) 2 ~ (1 -2nnb 3 ) c 2 k 2 + (vr/5) nb 5 c 2 k A + . . . . (2.16) 

Apparently, the result (12.161) for randomly orientated r = 2 defects agrees, to the order 
shown, with the previous result (12.141) for unoriented r = 1 defects. Some results for 
anisotropic defect distributions are given in App. [A] 
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For the r = 3 spacetime-foam model, the calculation is more complicated as the correction 
fields of the two "wormhole mouths" [spheres Sb and S' b in Eq. (12.61) ] affect each other. 
Therefore, we have to work directly with Taylor expansions in b/d. Giving only the leading 
order terms in b/d, the end result is 

^ r=3l (A;)) 2 ~(l-(20 7r/3)n6 3 ) c 2 k 2 + (2tt/9) nb 3 d 2 c 2 k" + . . . , (2.17) 

which holds for kb <C kd <C kl <C 1. For anisotropic defect distributions, some results are 
again given in App. [A] 

In closing, it is to be emphasized that any localized defect (weaving error) of space 
responds to an incoming electromagnetic plane wave by radiation fields corresponding to 
fictitious multipoles Only the position and relative strengths of these multipoles depend 
on the detailed structure of the defect. Together with the statistical distribution of the 
defects, these details then determine the precise numerical coefficients of the modified photon 
dispersion relation written as a power series in k (see Sec. II I II for further discussion). 



C. Scattering 

In this subsection, another aspect of electromagnetic-wave propagation is discussed, 
namely, the scattering of an incoming plane wave by r = 1 defects. Similar results are 
expected for r = 2 and r = 3' defects. 

As mentioned in Sec. Ill B[ the boundary conditions of a r = 1 defect correspond precisely 
to those of a perfectly conducting sphere. So the problem to consider is the scattering of 
an electromagnetic wave by a random distribution of identical perfectly conducting spheres 
with radii b and mean separation I, in the long-wavelength limit. More precisely, the relevant 
case has b ^ I ^ X, whereas ideal Rayleigh scattering (incoherent scattering by randomly 
distributed dipole scatterers) would have b <^ X <^ I . This means that all dipoles in a 
volume k~ 3 radiate coherently and their number, A^ co h ~ (k~ 3 )/(l 3 ), appears as an extra 
numerical factor in the absorption coefficient compared to standard Rayleigh scattering 
(see, e.g., Sec. 1-32-5 of Ref. 15 ] and Sees. 9.6 and 9.7 of Ref. (l6|). 



The relevant absorption coefficient (inverse scattering length) is then given by 

a£=; ] = l/L [ :=S ~ o-dip T 3 AU , (2.18) 

with (Tdip the cross section from the electric/magnetic dipole corresponding to an individual 
defect, l~ 3 the number density of such dipoles (i.e., defects), and N co ^ 3> 1 the coherence 
factor for the I <C A case discussed above. From the calculated polarizabilities of a r = 1 
defect, one has ~ k 4 b 6 neglecting factors of order unity. With N coh ~ (kl)~ 3 , the 
scattering length becomes 

~ k- 1 (i/6) a , (2.19) 

again up to factors of order unity. Expression (I2.19p suffices for our purpose but can, in 



principle, be calculated exactly, given the statistical distribution of defects [14|, [16 
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D. Proton dispersion relation 

In this last subsection, we obtain the dispersion relations from the Klein-Gordon and 
Dirac equations for the r = 1 spacetime-foam model. For the Klein-Gordon case, similar 
results are expected from the r = 2 and r = 3' models, but, for the Dirac case, the 
expectations are less clear and a full calculation seems to be required. 

For r = 1 defects and the long- wavelength approximation kb ^ 1 (i.e., considering the 
undisturbed harmonic fields to be spatially constant on the scale of the defect), the heuristics 
is as follows: 

(i) a scalar field obeying the Klein-Gordon equation does not require fictional sources to 
satisfy the boundary conditions at \x\ = b and, therefore, the dispersion relation is 
unchanged to leading order (there may, however, be other effects such as scattering 

H); 

(ii) a spinor field obeying the Dirac equation does require fictional sources but their con- 
tributions average to zero for many randomly positioned defects and the dispersion 
relation is unchanged to leading order. 

A detailed calculation (not reproduced here) gives, indeed, unchanged constant and 
quadratic terms in the dispersion relation of a real scalar, at least to leading order in kb. The 
Dirac calculation is somewhat more subtle and details are given in App. [B] The end result 
for the dispersion relation of a free Dirac particle, for definiteness taken to be a proton, is: 

(4 r=1] (A:)) 2 ~h~ 2 c 4 m 2 p + c 2 k 2 + ... , (2.20) 

with higher-order terms neglected and proton mass m p . These neglected higher-order terms 
in the proton dispersion relation would, for example, arise from additional factors k 2 b 2 and 
b 2 /l 2 , resulting in possible terms with the structure c 2 k 2 (m 2 c 2 /h 2 ) (b A /l 2 ) and c 2 fc 4 (b 4 /l 2 ). 

The combined photon and proton dispersion-relation results will be discussed further in 
the next section. 

III. DISPERSION RELATION RESULTS 
A. Coefficients and comments 

The different dispersion relations encountered up till now can be summarized as follows: 

(w? ) 2 = tr 2 c 4 m 2 s + (1 + K rs2 bl/ll) c 2 k 2 + K TS4 b 5 T /l 3 T c 2 k* + . . . , (3.1) 

for defect type r = 1, 2, 3 and particle species s = 7, 0, p corresponding to the Maxwell, 
Klein-Gordon, and Dirac equations, respectively. Four technical remarks are in order. First, 
the implicit assumption of fl3.ll) is that the same maximum limiting velocity c holds for all 
particles in the absence of defects (that is, for particles propagating in Minkowski spacetime). 
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TABLE I: Quadratic and quartic coefficients K in the photon dispersion relation (|3.ip . for s = 7, 
m-y = 0, and defect type r. 





•Kr72 




T = 1 


-2tt 


tt/5 


T = 2 


-2vr 


tt/5 


T = 3 


-20vr/3 


(2tt/9) d 2 /6 2 



Second, the photon mass vanishes in Maxwell theory, m 7 = 0, as long as gauge invariance 
holds. Third, only a few terms have been shown explicitly in ( 13. II) and, a priori, there may 
be many more (even up to order k A , as explained at the end of Sec. Ill Dj) . Fourth, a suffix r 
has been added to the length scales b and I of the models, since the length scale b of a r = 1 
defect, for example, is not the same quantity as the length scale b of a r = 2 defect. But, 
elsewhere in the text, this suffix is omitted, as long as it is clear which model is discussed. 

In the previous section, the photon coefficients K Tl2 and K T1 ± have been calculated for 
all three foam models (r = 1, 2, 3), but those of the scalar and proton dispersion relations 
only for the r = 1 model. The quadratic and quartic photon coefficients are given in Table [B 
The quadratic proton coefficient K lp % from the r = 1 foam model vanishes according to 
Eq. (I2.20p . as does the scalar coefficient K\^2 ■ Note that the present article considers only 
pointlike defects but that, in principle, there can also be linelike and planelike defects which 



give further terms in the modified dispersion relations [20[ . 

Let us close this subsection with three general comments. First, the modification of the 
quadratic coefficient of the photon dispersion relation, as given by Eq. (13.11) and Table [J 
can be of order unity (for b T somewhat less than l T ) and is not suppressed by powers of the 
quantum-electrodynamics coupling constant a or by additional inverse powers of the large 
energy scale A ~ hc/b T (which is already impossible for dimensional reasons, with m 7 = 
and a fixed density factor 1//^). This last observation agrees with a well-known result 



from quantum field theory; see, e.g., Refs. 2l|, |22|. Namely, if a symmetry (here, Lorentz 



invariance) of the quantum field theory considered is violated by the high-energy cutoff A (or 
by a more fundamental theory), then, without fine tuning, the low-energy effective theory 
may contain symmetry-violating terms which are not suppressed by inverse powers of the 
cutoff energy A. 

Second, the calculated dispersion relations (13.11) do not contain cubic terms in k, con- 
sistent with general arguments based on coordinate independence and rotational invariance 



23j. Furthermore, the photon dispersion relations found are the same for both polarization 
modes (i.e., absence of birefringence). For an anisotropic distribution of defects of type 
t = 2 or t = 3, however, the photon dispersion relations do show birefringence but still no 
cubic terms; see App. [S] 

Third, an important consequence of having different proton and photon dispersion rela- 
tions (13. ip is, as mentioned in the Introduction, the possibility of having so-called "vacuum 
Cherenkov radiation" 0,0]. A detailed study of this process in a somewhat different context 
(quantum electrodynamics with an additional Chern-Simons term in the photonic action) 



has been given in Refs. [2J, |25 
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B. General form 



The previous results on the dispersion relation (13. ip for the proton and photon can be 
combined and rewritten in the following general form: 

u 2 p = h~ 2 c\ m\ + c 2 p k 2 + 0(k 4 ) , (3.2a) 

uj 2 = (l + a 2 t/f)c 2 p k 2 + a i b 5 /tc 2 p k 4 + ..., (3.2b) 

for < kb, kl <C 1 and sign factors cf 2 , W4 G { — 1, 0, +1}. The velocity squared c 2 is defined 
as the coefficient of the quadratic term in the proton dispersion relation (I3.2a[) and the 
effective proton mass squared fn 2 is to be identified with the experimental value. 

With the results of Table fl] and Eq. (I2.20p . it is possible to get the explicit expressions 
for the effective length scales b and I in terms of the fundamental length scales b T and l T of 
the spacetime model considered. Specifically, the r = 1 spacetime-foam model has 

6~l(T 1/2 6i, 7~ (2tt)- 1/3 1(T 1/2 h, a 2 = -l, 0=4=1 ^ , (3.3) 

with radius b\ of the individual defects (identical empty spheres with antipodal points iden- 
tified) and mean separation l\ between the different defects. In other words, the effective 
and fundamental length scales of the r = 1 model are simply proportional to each other 
with coefficients of order unity. 

Similar results are expected for the r = 2 and r = 3' models, defined by Eqs. ( 12. 5ft 
and (12.71) . respectively. More generally, one could have a mixture of different defects (types 
t = 1, 2, 3', or other), with calculable parameters b, I, and cf 2i 4 in the photon dispersion 
relation (13.2bl) . 

For the purpose of this article, the most important result is that the effective length 
scales b and I of the photon dispersion relation ( I3.2bp are directly related to the fundamental 
length scales of the underlying spacetime model. This is a crucial improvement compared to 



a previous calculation of anomalous effects from a classical spacetime foam [26|, |27|, where 
the connection between effective and fundamental length scales could not be established 
rigorously. 

The parametrization ( 13. 2D for the isotropic case holds true in general and will be used 
in the following. Its length scales b and I will simply be called the average defect size and 
separation, respectively. Moreover, b/l and bj A ratios of order one will be allowed for, even 
though the calculations of Sec. Ill Bl leading, for example, to the identifications ( 13. 3p . are 
only valid under the technical assumptions 6/2 < 1 and b/X <C 1. In short, the proposal is 
to consider a modest generalization of our explicit results. 



IV. ASTROPHYSICS BOUNDS 



The discussion of this section closely follows the one of some previous articles [27|, |28|, 129 
which investigated modified dispersion relations from an entirely different (and less general 
[3o| ) origin. For completeness, we repeat the essential steps and give the original references. 
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Note also that, for simplicity, we focus on two particular "gold-plated" events but that other 
astrophysical input may very well improve the bounds obtained here. 



A. Time-dispersion bound 

The starting point for our first bound is the suggestion [lOj that the absence of time 
dispersion in a highly energetic burst of gamma-rays can be used to obtain bounds on 
modified dispersion relations (see, e.g., Refs. [3l|, |32] for subsequent papers and Ref. |33| for 
a review). 

From the photon dispersion relation (]3.2bl) . the relative change of the group velocity 
v g = du/dk between two different wave numbers k\ and ki is given by [34| : 



Ac 



ki,k 2 



v g (ki) - v g (k 2 



v g (ki 



(3/2) Ifcjf- k\\ b 5 /l 



(4.1) 



where Ac/c on the left-hand side is a convenient short-hand notation and where 6 and 7 on 
the right-hand side can be interpreted as, respectively, the average defect size and separation 
(see Sec. IIII Bl for further discussion). 

A flare of duration At from an astronomical source at distance D, with wave-number 



range k\ <C ki = k, 



7, max 



= E. 



7, max 



/(he), constrains the relative change of group velocity, 



Ac/c < cAt/D. Using (14. ip . this results in the following bound: 



T /7k 3/2 t 



(6/7) 



b < 



a/3/2 



he 



E. 



7, max 



1.2 x 10" 26 m 



cAt 
~D~ 



2.0 TeV 



1/2 



E, 



7, max 



At 
280 s 



1/2 



1.3 x 10 16 s 
D/c 



1/2 



(4.2) 



with values inserted for a TeV gamma-ray flare from the active galaxy Markarian 421 ob- 
served on May 15, 1996 at the Whipple Observatory 35, 36|. (The galaxy Mkn 421 has a 
redshift z ~ 0.031 and its distance has been taken as D = cz/H w 124 Mpc, for Hubble 
constant Hq w 75 km/s/Mpc.) The upcoming Gamma-ray Large Area Space Telescope 
(GLAST) may improve bound (14.21) by a factor of 10 4 , as discussed in App. A of Ref. 29 



B. Scattering bound 

It is also possible to obtain an upper bound on the ratio 6/7 by demanding the scattering 
length L to be larger than the source distance D or, better, larger than D/100 for an allowed 
reduction of the intensity by a factor F = exp(/) = exp(100). In other words, the chance 
for a gamma-ray to travel over a distance D would be essentially zero if L were less than 
-D/10 2 (see discussion below). 

The relevant expression for the scattering length L from r = 1 defects has been given in 
Sec. Ill Gl Here, we simply replace 6 and / in result (I2.19P by the general parameters 6 and 
7, again allowing for the case 6 ~ 7. Demanding L> D/f, then gives the announced bound: 
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{b/l) Z < v 7 / (K nmx Dy 1/2 

21 ( f \ 1/ V2.0TeV\ 1/2 / 3.8 x 10 24 m\ 1/2 . 

K 16x10 {W>) te) (—5—) ' (4 ' 3) 

using the same notation and numerical values as in Eq. (14. 2p . 

Strictly speaking, bound ( 14. 3 j) is useless if / is left unspecified. The problem is to decide, 
given a particular source, which intensity-reduction factor F = exp(/) is needed to be 
absolutely sure that its gamma-rays would not reach us if L were less than D/f. Practically 
speaking, we think that a factor F = exp(lOO) is already sufficient, but the reader can make 
up his or her own mind. More important for bound (14. 3 p to make sense is that one must be 
certain of the source of the observed gamma-rays and, thereby, of the distance D. For the 
particular TeV gamma-ray flare discussed here, the identification of the source as Mkn 421 
appears to be reasonably firm 35 . 



C. Cherenkov bounds 



A further set of constraints follows from the suggestion j^, [qJ that ultra-high-energy cosmic 
rays (UHECRs) can be used to search for possible Lorentz-noninvariance effects (or possible 
effects from a violation of the equivalence principle). The particular process considered here 
is vacuum Cherenkov radiation, which has already been mentioned in the last paragraph of 
Sec. [EDS 

From a highly energetic cosmic-ray event observed on October 15, 1991 by the Fly's Eye 
Air Shower Detector 37|, |38(, Gagnon and Moore 39j have obtained the following bounds 
on the quadratic and quartic coefficients of the modified photon dispersion relation |40l |: 



3 x 10 



-23 



< a 2 6 3 /f<3x 1(T 23 



(7 x l(T 39 m) : 



< a 4 6 5 // 3 



< 



(5 x l(T 38 m) : 



(4.4a) 
(4.4b) 



for length scales b, I and sign factors W2, as defined by Eqs. ( I3.2al) and ( I3.2bl) . For 
these bounds, the primary was assumed to be a proton with standard partonic distributions 
and energy ^uhecr ~ 3 x ]_qii GeV . Note that the limiting values of bounds (I4.4ap and 
(l44bj) scale approximately as (3 x 10 11 GeV/£" HECR ) n with n 



2 and n 



respectively. See Ref. 
heuristic discussion. 



39j for further details on these bounds and App. B of Ref. [29 



n b = 4, 
for a 



V. CONCLUSION 



The time- dispersion bound ( 14. 2ft on a particular combination of length scales from the 
modified photon dispersion relation (I3.2bl) is direct and, therefore, reliable. The same holds 
for the scattering bound fj4.3j) . as long as the allowed intensity-reduction factor is specified 
(see Sec. IIVBI for further discussion). The Cherenkov bounds (I4.4a[) and (I4.4b|) . however, 



are indirect in that they depend on further assumptions, e.g., interactions described by 
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quantum electrodynamics and standard-model partonic structure of the primary hadron. 
Still, the physics involved is well understood and, therefore, also these Cherenkov bounds 
can be considered to be quite reliable (ill ]. 

Turning to theoretical considerations, it is safe to say that there is no real understanding 



of what determines the large-scale topology of space [42]. With the advent of quantum 
theory, a similar lack of understanding applies to the small-scale structure of space 
Even so, assuming the relevance of a classical spacetime-foam model (see discussion below), 
Occam's razor suggests the model to have a single length scale, with average defect size b and 
average defect separation I of the same order (see Sec. IIIIBI for details on the interpretation 
of these length scales). Without natural explanation, it would be hard to understand why 
the defect gas would be extremely rarefied, b -C I. In the following discussion, we focus on 
the single-length-scale case but the alternative rarefied-g should be kept in mind. 

According to the time-dispersion bound (14.21) . a static classical spacetime foam with a 
single length scale (b ~ I) must have 



_ single-scale 
I 



< 1(T 26 m w he/ (2 x 10 10 GeV) , (5.1) 



which is a remarkable result compared to what can be achieved by particle-collider exper- 
iments on Earth. As mentioned in Sec. IIV A\ the experimental bound (15.11) may even be 
improved by a factor 10 4 in the near future, down to a value of the order of 10 -30 m. 

But the scattering and Cherenkov bounds (14.31) and f)4.4al) lead to a much stronger con- 
clusion: within the validity of the model, these independent bounds rule out a single-scale 
static classical spacetime foam altogether, 

b/1 <IQ- 7 . (5.2) 

The fact that a single-scale foam model is unacceptable holds even for values of b ~ I down 
to some 10~ 33 m « 10 2 x Zpi anc k (the precise definition of /pi an ck will be given shortly), at 
which len gth scale a classical spacetime may still have some relevance for describing physical 
processes 
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The unacceptability of a single-scale classical spacetime foam applies, strictly speaking, 
only to the particular type of models considered. But we do expect this conclusion to hold 
more generally (recall, in particular, the remarks of the last paragraph in Sec. Ill Bp . For 
example, also a time- dependent classical spacetime-foam structure with a single length scale 
appears to be ruled out 
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At distances of the order of the Planck length, Zpi anc k = yG h/c 3 ~ 1.6 x lCT 35 m f 
he/ (1.2 x 10 19 GeV), it is not clear what sense to make of a classical spacetime picture 45 
Still, at distances of order 10 2 x Zpi anc k, for example, one does expect a classical framework 
to emerge and, then, result (15.21) implies that the effective classical spacetime manifold is 



remarkably smooth [49|. If this conclusion is born out, it would suggest that either the 
Planck-length fluctuations of the quantum spacetime foam A a a, a 

are somehow made 

inoperative over larger distances or there is no quantum spacetime foam in the first place 



53 . 
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APPENDIX A: BIREFRINGENCE 

The individual defects of type r = 2 and r = 3 have a preferred direction given by the 
unit vector a in Figs. [2] and [31 respectively. An anisotropic distribution of defects may then 
lead to new effects compared to the case of isotropic distributions considered in the main 
text. This appendix presents some results for the photon dispersion relations from aligned 
r = 2 and r = 3 defects. 

Consider, first, a highly anisotropic distribution of r = 2 defects (empty spheres with 
points identified by reflection in an equatorial plane with normal vector a) having perfect 
alignment of the individual a vectors (henceforth, indicated by a caret, r = 2). Then, the 
two polarization modes (denoted © and 0) have dispersion relations: 



10. 



r=2] 



7* 



10. 



[r=2] 

79 



x2 (l + 67cnbh(kb) (k 2 -2k 2 )/k A ) c 2 k 2 

(ku,k ± )) - V ^-^ , (Ala) 

' (1 + 6nnb h{kb)/k 2 ) (l - 12nnb h{kb)/k 2 ) 

x2 (1 + 6nnb h(kb) (k 2 -2k 2 )/k 4 ) c 2 k 2 

- t \ , ,„, „ :\t J J \ f— - , (Alb) 



(1 + 6nnb h(kb)/k 2 ) (l - I2nnb h(kb)/k 2 ) 

with an auxiliary function h(x) = cosx — (sinx)/x = 0(x 2 ) for x — kb <C 1, parallel wave 
number k\\ = \k ■ a\, perpendicular wave number k± = \k — (k -a)a\, and defect number 
density n = l// 3 . For generic wave numbers with kn ^ k±, the two polarization modes have 
different phase velocities (v p hasc = ku/\k\ ) and there is birefringence. 

Consider, next, perfectly aligned r = 3 defects, that is, wormhole-like defects with two 
empty spheres identified by reflection in a central plane (normal vector a) and with all central 
planes parallel to each other (all vectors a aligned). In this case, we do not have a closed 
expression for the dispersion relations of the two polarization modes but rather Taylor series 
(the situation considered has kb kd kl 1). Only the results for two special wave 
numbers are given here. First, the photon dispersion relations for wave number parallel to 
the uniform orientation a of the defects (k± = 0) are equal for both polarization modes: 



(c4 T = 3] (%0)) =(4l^(%0)) ~(l-47rn6 3 ) c 2 k 2 + 2ir nb 3 d 2 c 2 k\ + ... . (A2) 

Second, the photon dispersion relations for wave number perpendicular to the defect ori- 
entation {k\\ = 0) are different for the two polarization modes (i.e., show birefringence): 

(J;f ] (0, fcj.)) 2 ~ (1 + 8tt nb 3 ) c 2 k 2 ± - (8tt/5) nb 5 c 2 k\ + . . . , (A3a) 
(J;f ] (0, k ± )Y ~ (1 - lQnnb 3 ) c 2 k\ + (16tt/5) nb 5 c 2 k\ + . . . , (A3b) 
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neglecting terms suppressed by a factor of b/d. 

The results of this appendix make clear that birefringence only occurs if there is some 
kind of "conspiracy" between individual asymmetric defects in the classical spacetime foam. 



APPENDIX B: DIRAC WAVE FUNCTION 



In this Appendix, we use the Dirac representation of the 7-matrices (for metric signature 

H and global Minkowski coordinates) and refer to Refs. 54], |55|, |56( for further details. 

For simplicity, we also set c = K — 1. The Dirac equation in Schrodinger form reads then 

id t ip(x,t) = (-1Q- V +m(3)ip(x,t), (Bl) 

with 4x4 matrices 



a 







7o7 



7o 





a 


c 





1 2 








-1 2 



(B2a) 



(B2b) 



in terms of the 2x2 unit matrix I2 and the 2x2 Pauli matrices a. 

In the presence of a single r = 1 defect (sphere with radius b centered at x = and 
antipodal identification), we impose the following boundary condition on the Dirac spinor: 

[T=l] 



if>(—x, t) = 1 x ■ a i/)(x, t) 



\x\=b 



(B3) 



for unit vector x = x/\x\. [There can be an additional phase factor rj G C (\r]\ = 1) on the 
right-hand side of (1B3j) . which may in principle depend on the direction, 77 = f](x).} The 
physical motivation of boundary condition (1B3[) is that the Dirac particle moves appropri- 
ately near the defect at |af| = b. Recall that, in the first-quantized theory considered (cf. 
Ref. |55j|). the free particle is described by a wave packet which can in principle have an 
arbitrarily small extension, for example, much less than b. The boundary condition (IB3j) 
then makes the probability 4-current j M = ip^^ip = ijr (I4, a) if) well behaved near the defect 
at \x\ = b: probability density j° equal at antipodal points, normal component of j going 
through, and tangential components of j changing direction (cf. Fig. [1] and the discussion 
in Sec. EH). 



The case of primary interest to us has spin-| particles of very high energy compared to 
the rest mass m but wavelength still much larger than the individual defect size b and the 
mean defect separation 



m<A;<l/Z<l/&. 
An appropriate initial solution of the Dirac equation over M 4 is given by 

/l\ 


1 

V / 



(B4) 



t)in ~ exp (1 kz — 1 ut) 



(B5) 
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which corresponds to a positive-energy plane wave propagating in the z = x 3 direction. This 
wave function, however, does not satisfy the defect boundary condition (IB3I) for the 3-space 
(12. 2p with a single defect centered at the point x — 0. 

Make now the following monopole-like Ansatz for the required correction: 



ip{x,t) ( 



exp (—1 cut) g(r/b)\x-a 



S3 
\S,J 



(B6) 



with radial coordinate r = \x\, normalization g{l) = 1, and complex constants s n . (This par- 
ticular Ansatz is motivated by the structure of the Green's function for the Dirac operator; 
cf. Sec. 34 of Ref. 1541 .) The total wave function, 



(B7) 



■ifj(x, t) = ip(x, t) in + ip(x, t) 



must then satisfy the defect boundary condition (IB3I) at r = b, in the limit kb — * 0. The 
appropriate constant spinor (si, S2, S3, S4) in (1B6I) is readily found. Also, the function g in 
the "near zone" (r C A) must be given by g{r/b) = b 2 /r 2 , in order to satisfy the Dirac 
equation neglecting terms of order mb, kb, and m/k. 

All in all, we have for the corrected wave function from a single defect centered at x = x±: 



ij) (x, t) ~ exp (1 kz — 1 ut ) 



(l\ 

1 

V / 



1 exp ( — 1 C(jt) 



/ cos 61 \ 
sin#! exp(i0!) 
cos 6*i 
V sin 6*i exp(i0i) J 



(B8) 



where (ri,9i,<pi) are standard spherical coordinates with respect to the defect center x,\ and 
the z axis from the global Minkowski coordinate system, having r\ = \x—xi \ > b, 9\ G [0, ir), 
and 0i G [0,2n). 

Next, sum over the contributions of identical defects with centers x = Xj, for j = 
1, . . . , N. The resulting Dirac wave function at a point x between the defects is given by 



ip(x, t) ~ exp (1 kz — 1 ut) 




1 

V / 



1 exp 



-10J 



N 

i=i 



/ cos 9j \ 
sm9j exp(i0j) 
cos 9j 
y sm9j exp(i0j) J 



(B9) 



where 9j and <f)j are polar and azimuthal angles with respect to the defect center Xj and 
the z axis. With many randomly positioned defects present (N ^> 1), the entries of the 
second spinor on the right-hand side of (1B9j) average to zero and only the initial Dirac 
spinor remains. 

For the electromagnetic case discussed in Sec. Ill Bl the fictional electric/magnetic dipoles 
(radial dependence oc 1/r 3 ) are aligned by the linearly polarized initial electric/magnetic 
fields (12 .9p and there remain correction fields after averaging, which produce a modification 
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of the photon dispersion relation. As mentioned above, the corrective wave function required 
for the Dirac spinor is monopole-like (radial dependence oc 1/r 2 ) and averages to zero. This 
different behavior is a manifestation of the fundamental difference between vector and spinor 
fields. 

The final result is that the dispersion relation of a high-energy Dirac particle is unchanged, 
uj 2 ~ k 2 , at least up to leading order in m/k, kl, and b/l. For an initial spinor wave function 
at rest (k = and & C ! < 1/m), a similar calculation gives uo 2 ~ m 2 . (Note that this is 
an independent calculation as Lorentz invariance does not hold.) Combined, the dispersion 
relation of a free Dirac particle (for definiteness, taken to be a proton) reads 

(^ =1 l) 2 ~m 2 + fc 2 + ... , (BIO) 

to leading order in ml, kl, and b/l. 

To summarize, we have found in this appendix that the classical spacetime-foam model 
considered affects the quadratic coefficient of the proton dispersion relation differently than 
the one of the photon dispersion relation. However, the calculation performed here was only 
in the context of the first-quantized theory and a proper second-quantized calculation is left 
for the future. 
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